Contact complete integrability 

Boris Khesin* and Serge Tabachnikov^ 
October 1, 2009 



Abstract 

Complete integrability in a symplectic setting means the existence 
of a Lagrangian foliation leaf-wise preserved by the dynamics. In 
the paper we describe complete integrability in a contact set-up as 
a more subtle structure: a flag of two foliations, Legendrian and co- 
Legendrian, and a holonomy-invariant transverse measure of the for- 
mer in the latter. This turns out to be equivalent to the existence of a 
canonical M c< M™ -1 structure on the leaves of the co-Legendrian foli- 
ation. Further, the above structure implies the existence of n contact 
fields preserving a special contact 1-form, thus providing the geometric 
framework and establishing equivalence with previously known defini- 
tions of contact integrability. We also show that contact completely 
integrable systems are solvable in quadratures. 

We present an example of contact complete integrability: the bil- 
liard system inside an ellipsoid in pseudo-Euclidean space, restricted 
to the space of oriented null geodesies. We describe a surprising accel- 
eration mechanism for closed light-like billiard trajectories. 
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1 Introduction 

Our first motivation is the following V. Arnold's problem No 1995-12 in [I]: 

Transfer the theory of completely integrable Hamiltonian sys- 
tems from symplectic geometry to contact geometry (where, e.g., 
the Lagrangian invariant manifolds with their natural affine struc- 
tures determined by Lagrangian fibrations must be substituted 
by Legendrian invariant manifolds with their natural projective 
structures determined by Legendrian fibrations). Carry over the 
Liouville theorem to this context and find applications to the 
infinite-dimensional case (where the equations of characteristics 
are partial differential). 

The classical set-up for the Arnold-Liouville theorem is a symplectic 
manifold (M 2n ,uj) (for example, the phase space of a mechanical system) 
and a discrete- or continuous-time symplectic dynamical system on it, that 
is, a symplectomorphism T : M — > M or a symplectic vector field v on M, 
respectively. (Here and elsewhere we refer to (2j for a succinct exposition 
of the basic facts of symplectic and contact geometry.) Recall that a La- 
grangian manifold F n C M 2n is a half-dimensional submanifold such that 
the restriction of u to F vanishes. A symplectic dynamical system is called 
completely integrable if M is endowed with a Lagrangian foliation T whose 
leaves are invariant under the dynamics. 

A fundamental geometrical fact, underlying the Arnold-Liouville theo- 
rem, is that the leaves of a Lagrangian foliation carry a canonical affine 
structure. Choose n functionally independent "integrals" (functions, con- 
stant on the leaves of J-) and consider their symplectic gradients. One 
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obtains n commuting vector fields, tangent to the leaves of T and providing 
a field of frames along the leaves. A different choice of integrals results in 
applying linear transformations, constant along each leaf, to these frames. 
Thus each leaf has a flat structure. 

The map T, or the vector field v, preserves the symplectic structure and 
the foliation leaf-wise, and hence preserves the afhne structure on the leaves. 
It follows that T is a parallel translation, and v is a constant vector field, 
on each leaf of the Lagrangian foliation. 

This has strong dynamical consequences. If a point is periodic then so are 
all the points on the same leaf of and with the same period. This implies 
Poncelet-style theorems (see, e.g., |14] for a recent application). If a leaf is 
compact, it must be a torus, and the dynamics is a quasi-periodic motion 
on the torus. Another useful consequence: if two symplectic maps share an 
invariant Lagrangian foliation then they commute (because so do parallel 
translations). We refer to [22] for more detail on complete integrability in 
the discrete-time case. 

To summarize, the definition of a completely integrable dynamical sys- 
tem consists of two parts: a certain geometrical structure on a symplectic 
manifold M, namely, a Lagrangian foliation J 7 , and a discrete- or continuous- 
time symplectic dynamical system, preserving this structure. It is natural 
to call the first part, the pair (M, J 7 ), a completely integrable symplectic 
manifold. 

Contact manifolds are odd-dimensional relatives of symplectic manifolds. 
Let (M 2n-1 ,£) be a contact manifold with a contact distribution £. Recall 
that a Legendrian submanifold F n ~ l C M 2n ~ l is an integral manifold of £ of 
the maximal possible dimension n — 1. The leaves of a Legendrian foliation 
carry a canonical projective structure: this is a contact counterpart to the 
above-mentioned theorem about Lagrangian foliations (we shall dwell on 
this projective structure in Sections 2.2). The problem is to extend the 
notion of complete integrability to contact manifolds. 

Note that the simplest particular case of a contact manifold is a 1- 
dimensional manifold, M 1 or S 1 , with the trivial contact structure. A natural 
definition of integrability in dimension 1 (with discrete- or continuous-time) 
is the existence of a non-vanishing invariant differential 1-form; for a diffeo- 
morphism of S , this implies that the map is conjugated to a rotation. 

Contact complete integrability was studied before: see [5J El [131 US] 
and also related papers on Legendrian foliations |19[ [T5] . For example, 
according to [5 , a completely integrable vector field on a contact manifold 
M 2n ~ 1 is the Reeb field of a contact 1-form, for which the space of first 
integrals determines a fibration with n-dimensional fibers defined locally by 
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the action of a torus T n of contact transformations, see Section 3.5 for a 
brief survey of earlier work. 

The main goal in this paper is to give a definition that extends the earlier 
ones and that works equally well in the continuous- and discrete-time cases. 
Our second motivation was to place the recently studied examples j9j [TO] 
into the general context of contact complete integrability. 

These examples include the geodesic flow on an ellipsoid and the billiard 
map inside an ellipsoid in pseudo-Euclidean space. In pseudo-Euclidean 
setting, one has a trichotomy for an oriented line: it may be space-like, 
time-like, and light- like (or null), that is, having positive, negative, or zero 
energy. The manifolds of oriented non-parameterized space- and time-like 
lines carry canonical symplectic structures, just like in the Euclidean case, 
but the space of null lines has a canonical contact structure; see [TO] and 
Section |U 

Let S be a smooth closed hypersurface in a pseudo-Euclidean space. 
The billiard system inside S can be considered as a map on the space of 
oriented lines taking the incoming billiard trajectory to the outgoing one. 
The law of reflection is determined by the energy and momentum conser- 
vation, therefore the type of a line (space-, time-, or light-like) does not 
change. Restricted to space- and time-like lines, the billiard transformation 
is a symplectic map, but its restriction to light-like lines is a contact map. 

If S is an ellipsoid, the respective billiard transformation is integrable, 
in the following sense. An ellipsoid in n + 1-dimensional pseudo-Euclidean 
space determines a pseudo-confocal family of quadrics, see [TO] and Sec- 
tion [4] A space- or time-like billiard trajectory remains tangent to n fixed 
pseudo-confocal quadrics. This gives n integrals of the billiard map on the 
2n-dimensional symplectic spaces of oriented space- or time-like lines. These 
integrals Poisson commute and hence define an invariant Lagrangian folia- 
tion. This is just like the Euclidean case, see, e.g., [20| [21] . 

However, we lose one integral on the space of null lines: a light-like 
billiard trajectory remains tangent to n — 1 fixed pseudo-confocal quadrics. 
This gives n—1 integrals on the In — 1-dimensional contact space of oriented 
light-like lines and hence a foliation T n . It turns out that the distribution 
given by the intersection of the tangent spaces to the leaves of T with 
the contact hyperplanes is also integrable, and one obtains a Legendrian 
foliation Q n ~ l whose leaves foliate the leaves of T. Furthermore, the billiard 
transformation has an invariant contact form - morally, another integral, 
since all contact forms for a given contact structure differ by multiplication 
by a non-vanishing function - and this additional integral commutes, in an 
appropriate sense, with the other n — 1 integrals, see Section [4] 
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The above motivates the following general definition of contact integra- 
bility. Let (M 2n_1 ,^) be a contact manifold with contact distribution £. A 
foliation T n is called co-Legendrian if it is transverse to £ and the distribu- 
tion TJ^n^ is integrable. Let Q n ~ l be the respective Legendrian foliation. 
We have a flag of foliations (J-, Q). In Section [2] we show that the canonical 
projective structure on the leaves of Q reduces to an afhne structure. 

If a contact dynamical system preserves a co-Legendrian foliation leaf- 
wise then it sends the leaves of the respective Legendrian foliation Q to each 
other, preserving the affine structures therein. Thus the dynamics reduces 
to 1-dimensional one on the space of leaves of Q within a leaf of T '. For 
this dynamics to be integrable, one needs an invariant 1-form on this 1- 
dimensional space of leaves. 

Definition 1 A completely integrable contact manifold M is a co-Legendrian 
foliation fonM such that, for each leaf F of J 7 , the respective codimension 
one foliation Q on F has a holonomy invariant transverse smooth measure. 
A discrete- or continuous-time contact completely integrable system on M 
is a contactomorphism, or a contact vector field, that preserves J- leaf-wise 
and preserves the above transverse measure of the foliation Q. 



We show in Section 3.1 that the leaves of a co-Legendrian foliation on a 
completely integrable contact manifold have a canonical M. ix M n_ ^structure. 
This has strong dynamical implications, similarly to the flat Restructure 
on the leaves of a Lagrangian foliation of a symplectic manifold. 

An example of a completely integrable contact manifold M is analyzed 
in Sections |3.2| and |3.3| M has a contact form whose Reeb field is tangent 
to the co-Legendrian foliation. We show that then the contact form defines 
a holonomy invariant transverse smooth measure on the space of leaves of 
the Legendrian foliation within a leaf of the co-Legendrian one. We show 



in Section 3.4 that the familiar example of a completely integrable geodesic 
flow on a Riemannian manifold fits into this framework. 

In Section [4] we show that the billiard ball map inside an ellipsoid in 
pseudo-Euclidean space, restricted to oriented light-like lines, is a completely 
integrable contact map. We do this by constructing an invariant contact 
form on the contact space of oriented null lines whose Reeb field is tangent 
to a co-Legendrian foliation. 



2 Geometry of co-Legendrian foliations 

In this section we study the geometry of co-Legendrian foliations. 
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2.1 Example of a co-Legendrian foliation 



Example 2.1 Let M be a contact manifold with a Legendrian foliation Q, 
and let <fit be a 1-parameter group of contactomorphisms preserving this 
foliation. Assume that the vector field corresponding to 4>t is transverse to 
the contact distribution. Then, acting by <pt on Q, yields a co-Legendrian 
foliation J 7 , that is, the leaves of T are the orbits of the leaves of Q under 
the flow <f>t. 

In fact, this example is universal, as the following lemma shows. 
Lemma 2.2 Every co-Legendrian foliation is locally contactomorphic to the 



one in Example 2.1 



Proof. Recall that a contact element on a smooth manifold N is a hyper- 
plane in a tangent space to N. Since a contact element is the kernel of a 
covector, uniquely determined up to a non-zero factor, the space of contact 
elements is PT*N, the projectivization of the cotangent bundle. The space 
of contact elements has a canonical contact structure given by the "skat- 
ing condition": the velocity of the foot point of a contact element lies in 
this contact element. One has the fibration p : PT*N — > N whose fibers 
are Legendrian submanifolds; these fibers consist of contact elements with 
a fixed foot point. Every Legendrian foliation is locally contactomorphic to 
this one. 

Suppose that M is a contact manifold, J 7 is a co-Legendrian foliation 
and Q the respective Legendrian foliation. We may assume that, locally, 
M = PT*N and Q is the fibration p : PT*N N. Then T projects to 
a 1-dimensional foliation C in N. In other words, a leaf of J- consists of 
contact elements whose foot points lie on a leaf of C. 

Consider a 1-parameter group of diffeomorphisms of N whose trajecto- 
ries are the leaves of C Diffeomorphisms of N naturally act on contact 
elements on N, so we obtain a 1-parameter group of contactomorphisms 
of M preserving the foliation Q. If we restrict to the open set of contact 
elements on N that are not tangent to C then T is obtained from Q as in 
Example |2.1| □ 



2.2 Symplectic interpretation of co-Legendrian foliations 

Let (M 2n-1 ,£) be a contact manifold. Recall the notion of symplectization 
(or the symplectic cone). Let P 2n C T*M be the total space of the principle 
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M*-bundle ir : P —* M whose fibers consist of non-zero covectors (x,p) that 
vanish on the contact element in M at their respective foot points 
x. The symplectization P has a canonical 1-form A, the restriction of the 
Liouville 1-form in T*M, and the 2-form co = dX is a symplectic structure 
on P. One has the multiplicative M*-action on P; let 2? be the respective 
vector field, called the Euler field. The following identities hold: 

i E uJ = \, X(E) = 0, L E X = X. (1) 

For example, the symplectization of the space of contact elements PT*N 
is the cotangent bundle T*N with the zero section removed. The M*-action 
is the fiber- wise scaling of covectors, and the 1-form A is the Liouville form 
in T*N. 

The operation of symplectization relates the contact geometry of M 
to the homogeneous symplectic geometry of P. Specifically, contactomor- 
phisms of M are the symplectomorphisms of P that commute with the In- 
action; the preimage of a Legendrian submanifold in M is an M*-invariant 
(conical) Lagrangian submanifold in P; the preimage of a Legendrian folia- 
tion in M is an M*-invariant Lagrangian foliation in P, etc. 

Let T be a co-Legendrian foliation on M and Q the respective Legendrian 
foliation. Set: T = 7r _1 (J r ), Q = tt~ 1 {G). In the next lemma, we interpret 
co-Legendrian foliations in symplectic terms. 

— n+l 

Lemma 2.3 T is a co-isotropic foliation in P. Its symplectic orthogo- 
nal complement 7i. n is an isotropic foliation transverse to E, and Q n is 
spanned by E and H" (see figureVlv. Conversely, given a co-isotropic fo- 

— n+l 

Ration T in P, tangent to the Euler field E and transverse to ker A, the 
projection of J- to M is a co-Legendrian foliation therein. 

Proof. Let fa : M — ► R, i = 1, ... ,n — 1, be locally defined functions 
whose common level surfaces are the leaves of and let fa = 7r*(/»). Then 
the homogeneous functions fa : P — > K of degree zero define the foliation T . 

First, we show that the symplectic orthogonal complement to the tangent 
space TT is spanned by the Hamiltonian vector fields sgrad fa. Indeed, 
consider a vector v G TT. Then ^(sgrad/j, v) = dfa(v) = 0, since fa is 
constant on the leaves of T . 

Next, we show that the distribution spanned by the Hamiltonian vector 
fields sgrad fa is integrable. Indeed, this distribution is isotropic, hence 
w(sgrad/j, sgrad /,) = = {fa, fj}. It follows that [sgrad fa, sgrad fj] = 0, so 
7i is a foliation. 
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Figure 1: A leaf F of the foliation F in P, foliated by Q and by H projects 
to a leaf F of the foliation F in M, foliated by Q 

One has: Q = FC\ ker A. We claim that Q is spanned by E and sgrad/j. 
Indeed, E and sgrad/j are tangent to F. One has: X(E) = and 

A(sgrad/i) = u(E, sgrad/j) = -dfi(E) = 0, 

since fa is homogeneous of degree zero with respect to the Euler field. Thus 
E and all sgrad/j are tangent to Q. 

Let us check that E is transverse to H. If not, then, at some point, 
E = sgrad/ for a function / : M — ► R that is constant on the leaves of F . 
Then at that point A = Ieoj = i sgra d/ w = df. This is a contradiction since 
the foliation T is transverse to the contact structure, and hence A does not 
vanish on the tangent spaces to its leaves. 

— n+l 

Finally, we claim that if T is a co-isotropic foliation in P, tangent to 
the Euler field E, then the projection it : P — > M takes F to a co-Legendrian 
foliation. Indeed, the foliation F is invariant under the Euler field since E 
is tangent to it. Thus F is conical. Then the distribution TF n ker A is a 
conical Lagrangian foliation that projects to a Legendrian foliation in M. □ 

Thus a co-Legendrian foliation on a contact manifold M 2ri_1 is the same 
as a co-isotropic n + 1- dimensional foliation on its symplectization P 2n given 
by n — 1 Poisson commuting homogeneous functions of degree zero. 
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2.3 Flat structure on the leaves of Q 

As we mentioned in Introduction, the leaves of a Legendrian foliation carry 
a canonical projective structure. Let us recall this construction. 

Let (M, £) be a contact manifold and Q a Legendrian foliation. As before, 
we may assume that M = PT*N and Q is the fibration p : PT*N — ► AT. 
Let x & N and G x = Then dp takes the contact hyperplanes 

along the leaf G x to hyperplanes in the tangent space V := T X N. The 
set of all such hyperplanes is P(V) = MP" -1 , and we obtain a mapping 
ip : G x — > KP n_1 . Due to complete non-integrability of the contact structure, 
if is a local diffeomorphism. Thus G x has a projective structure. 

Now let (M,£) be a contact manifold, and J 7 and £/ be co-Legendrian 
and the respective Legendrian foliations. 

Lemma 2.4 T/ie projective structure on the leaves of Q has a reduction to 
an affine structure. 

Proof. In the notation of the preceding paragraphs, the tangent spaces to a 
leaf of T are taken by dp to a line £ C V. The set of hyperplanes in V passing 
through £ is a projective hyperplane MP™ -2 C P(V) = MP™ -1 , and the 
image of tp does not intersect this projective hyperplane. The complement 
j^pn-i _ - g an a fg ne S p ace . Thus we have a local diffeomorphism 

<p : G x — > A n_1 , whence an affine structure on the leaves of Q. □ 

As usual, the existence of an affine structure imposes restrictions on the 
topology of the leaves. For example, a compact leaf of Q is a torus. 

Remark 2.5 Alternatively, one can define an affine structure on the leaves 
of Q as follows. Recall that the isotropic foliation TC is generated by the 
vector fields sgrad/j, where the functions fi are homogeneous of degree zero. 
The commuting vector fields sgrad/j define an affine structure on the leaves 
of H. 

Since the functions fi are homogeneous of degree zero, [E, sgrad/j] = 
— sgrad/j. Therefore the M*-action preserves the foliation 7i, sending leaves 
to leaves, and these maps preserve the affine structure on the leaves. The 
projection it : P — > M diffeomorphically maps the leaves of TC to the leaves 
of the Legendrian foliation Q endowing the latter with an affine structure. 

2.4 Weakly integrable contact systems 

Definition 2 A discrete- or continuos-time contact weakly integrable system 



9 



is a contact dynamical system on a contact manifold that has a leaf-wise 
invariant co-Legendrian foliation. 



Such a system reduces to a 1-dimensional one. The leaves of the Leg- 
endrian foliation Q within a leaf F of the co-Legendrian foliation J- are 
mapped to each other by parallel translations in their respective affine co- 
ordinates, but the motion on the 1-dimensional space of leaves F/Q may 
be arbitrary. As Lemma 2.2 shows, any diffeomorphism of N, preserving 
the one- dimensional foliation C leaf-wise, lifts to a weakly integrable con- 
tactomorphism of (an open subspace of) the space of contact elements of 
N. 

If M is 1-dimensional, the co-Legendrian foliation consists of one leaf, 
M itself, and the definition imposes no constraints on the dynamics. 



3 Completely integrable contact manifolds 

In this section we study the geometry of completely integrable contact man- 
ifolds and completely integrable contact dynamical systems. 



3.1 Semi-direct product structure 

Let G be a subgroup of the group of diffeomorphisms of R n . A G-structure 
on an n-dimensional manifold is (an equivalent class of) an atlas whose 
transition functions belong to G. In these terms, the leaves of a Lagrangian 
foliation have an Restructure, where R n is the group of parallel translations 
of n-dimensional affine space. 

Let M 2n_1 be a completely integrable contact manifold with the flag 
of co-Legendrian and Legendrian foliations (J- T \G n ~ 1 ). Let R x W l ~ l be a 
semi-direct product of R and R n_1 : 

K"- 1 -4R« R"- 1 R 0. 

Example 3.1 Given a number A £ R and a vector b 6 R™ _1 consider 
affine maps v i— > e x v + b of the space R n_1 9 Then the set of such 
pairs (A, b) forms a Lie group with respect to natural composition of the 
affine maps. Similarly, one can define the Lie group by composing affine 
transformations v \— > e x Pv + b for a projector P : R n_1 — > R 71 ^ 1 (with 
P 2 = P). These Lie groups give examples of semi-direct products R >< R n_1 . 
The first example corresponds to the case of P = id, while the direct product 
group R n = R x R n_1 corresponds to P = 0. 
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Lemma 3.2 The leaves of the co-Legendrian foliation J 7 of a completely- 
integrable contact manifold have a canonical R k M n_1 -structure. 



Proof. Let F be a leaf of J- . One has an exact sequence of vector bundles: 
— > TQ — > TF —> NQ = TF/rg -> 



(iVC? is the normal bundle of the foliation Q). By Lemma 2.4 the leaves of 
Q have an M n_1 -structure. The transverse invariant measure of the foliation 
Q in F fixes a trivialization of NQ. The two combined yield a M k IR n_1 - 
structure on F. □ 

Recall that a completely integrable system on a symplectic manifold M 2n 
can be defined by a local IR n -action which preserves the symplectic structure 
and is generically free. 

Similarly, a completely integrable system on a contact manifold can be 
defined by a R >< M n ~ 1 -action, where the abelian subgroup IR n_1 acts locally 
free along the contact planes. Namely, consider a contact manifold M 2n_1 



with a contact distribution £. One can see that Lemma 3.2 is equivalent to 
the following 

Lemma 3.3 The existence of a co-Legendrian foliation with an invariant 
measure on a contact manifold (M,^) is equivalent to the existence of a local 
M k W l ~ 1 -action on M such that the W 1-1 -orbits of the abelian subgroup are 
tangent to the distribution £. 

By construction, the R >< M n-1 -orbits define leaves of the co-Legendrian 
foliation J 7 , while the orbits of the abelian subgroup provide the Legendrian 
foliation Q. The M-action in the semi-direct product gives the holonomy- 
invariant transverse measure. Conversely, the existence of the R ix R n_1 - 
structure for T implies the existence of a local M k M n_1 -action in each leaf 
of whose M ra_1 -orbits are leaves of the Legendrian foliation Q. 



3.2 A special contact form 

It turns out that the existence of a local M >< R n ^ 1 -action on a contact 
manifold M implies the existence of a special contact form whose Reeb field 
is tangent to the group orbits. Recall that the Reeb vector field v of a contact 
form A spans the kernel of dX and is normalized by the condition A(t> ) = 1. 

Lemma 3.4 There is a contact form A on M whose Reeb field is tangent 
to the E k W 1 ' 1 -orbits. 
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Proof. Let V £ (IRixIR™ -1 ) be a generic element of the Lie algebra and U € 
M n_1 an element of the abelian subalgebra. Let v and u be the corresponding 
vector fields on M. Note that in the semi-direct product R K W l ~ l the 
elements U and V (and hence the fields u and u) satisfy the relation [v, u] = 
au for some a £ R. Also note that by the definition of the action, u is 
tangent to the distribution £, while v is transversal to £ at a generic point. 

Now let Ao be a contact 1-form defining the contact structure £. Define 
the 1-form A by normalizing Ao as follows: A = Ao/Ao(i>), so that X(v) = lQ 

We see that, on the one hand, 

A = i a uX = . 

On the other hand, 

*[«,it] ^ = ivL u \ — L u i v \ = i v i u d\ + i v di u \ — L u i v \ = i v i u d\ , 

where we used that i^A = 1 and i u X = 0. The equality i v i u d\ = implies 
that the kernel of the 2-form dX is tangent to the orbits of the M ix M n_1 - 
action. Indeed, this equality shows that the projection of v to planes of £ 
along this kernel must be tangent to the Legendrian M n_1 -orbits. Thus the 
Reeb field for the 1-form A belongs to the M k M ra_1 -orbits. □ 

The following lemma shows that the converse statement also holds: the 
existence of such a special contact 1-form is equivalent to the existence of 
co-Legendrian foliation with a transverse measure. 

Let (M, £) be a contact manifold with a contact form A, and let T and Q 
be a co-Legendrian and the respective Legendrian foliations. Does a contact 
form A determine a holonomy invariant transverse smooth measure of the 
foliation Q within the leaves of T1 The next lemma also gives a sufficient 
condition. 

Lemma 3.5 Assume that the Reeb field of X is tangent to the co-Legendrian 
foliation T ' ■ Let F be a leaf of T . Then the 1-form X determines a holonomy 
invariant transverse smooth measure of the foliation Q on the manifold F. 

Proof. We need to check that the restriction of A to F is a basic differential 
form with respect to the foliation Q; this means that for every vector field 
u, tangent to Q, one has: i u X = L U X = 0. If A is basic then it descends on 
the (locally defined) space of leaves and defines a 1-form on this space. 

1 This construction of the invariant 1-form is similar to the one for the action of an 
abelian group discussed in [6]. 
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We have A(u) = 0, since Q is Legendrian. Then, by Cartan's formula, 
L U X = i u d\, and we want to show that i u dX = 0. The tangent space TF is 
spanned by TQ and v, the Reeb field. If w € TQ then dX(u, w) = since 
T£? is a Lagrangian subspace of the symplectic space £ = ker A with the 
symplectic structure dA. On the other hand, d\(u, v ) = since u 6 ker A. 
Thus i u d\ = 0, and we are done. □ 

The above two lemmas give a necessary and sufficient condition of con- 
tact integrability in terms of a special 1-form. 

Theorem 3.6 The existence of a co-Legendrian foliation T with an in- 
variant transverse smooth measure is equivalent to the existence of a co- 
Legendrian foliation with a special contact 1-form whose Reeb field is tangent 
to the foliation. 

3.3 Commuting fields and invariant contact forms 

The existence of a local R k M n_1 -action also implies the existence of an 
appropriate local M n -action (with the same orbits), preserving this form: 
one can define n commuting vector fields which leave the contact form A 
invariant and span the same foliation T . Note however, that although the 
M K M n_1 - and R n -orbits coincide, the orbits of the ^subgroups in these 
two groups are different: any nonzero vector field preserving a contact form 
cannot be tangent to a contact distribution. 

Recall that if a contact form A is chosen on a contact manifold (M, £) 
then one can assign a contact vector field Xf to a smooth function /: a 
contact form determines a section of the symplectization P, and this makes 
it possible to extend / to P as a homogeneous of degree one function; the 
Hamiltonian vector field of this extended function projects to a contact vec- 
tor field Xf on M. The correspondence between the functions and contact 
vector fields is described by the formula X(Xf) = f. In particular, for / = 1 
one has X\ = v, the Reeb field. Note also the formula: Lx f X = df(v)X. 

Further, one defines the Jacobi bracket on smooth functions: [/, g] = 
X([Xf, X g ]). This operation satisfies the Jacobi identity, but not the Leibniz 
one. The correspondence / i— ► Xf is a Lie algebra homomorphism. One has 
the identity: 

[f,g] = dX(X f ,X g ) + fdg(v) -gdf(v). (2) 

One also has a projection TM — ► £ along the direction of the Reeb field v. 
Denote by u the "horizontal" part of u G TM, that is, its projection to the 
contact hyperplane. Then one has: Xf = fv + Xf. 
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As in Lemma 3.5 assume that the Reeb field v is tangent to a co- 
Legendrian foliation T on a contact manifold M 2 ™" 1 with a contact form 
A. Let fi : M — ► K, i = 1, . . . , n — 1, be locally defined functions whose 
common level surfaces are the leaves of and let tij = X* v i = 1, . . . , n — 1. 

Lemma 3.7 TTie vector fields v, u±, . . . , it n — l pairwise commute and span 
the foliation T . 



Proof. Let / be an "integral" of the foliation T , that is, a function con- 
stant on the leaves. The formulas X(Xf) = f and Lx f X = df(v)X, along 
with the Cartan formula, imply that ix f dX = df(v)X — df. It follows that, 
for every test vector w G TQ, one has dX(Xf,w) = 0: indeed, X(w) = 
since w £ £, and df(w) = since w G TT . It follows that Xf lies in the 
symplectic orthogonal complement to TQ in £. Since £/ is Legendrian, Xf is 
tangent to ^. Therefore Xf is tangent to T . 

Next, we claim that [fi, fj] = 0. Indeed, since v is tangent to J-, one 
has dfi(v) = for all i. It follows from ^ that = dX(ui,uj). Since 

w G keriiA, the latter is equal to dX(ui,Uj), and this is zero since all lie 
in the Legendrian space TQ C £. 

Likewise, = dX(v,Ui) = since u G kerdA. It follows that the 

vector fields v, u±, . . . , u n -\ pairwise commute, as claimed. □ 

Suppose that a system of differential equations is given. To solve the 
system in quadratures means to obtain its solution by a finite number of 
"algebraic" operations (including inversion of functions) and "quadratures" , 
integration of known functions, see, e.g., [3]. 

Theorem 3.8 A continuous-time completely integrable contact system is 
solvable in quadratures. 



Proof. If the contact structure and n — 1 first integrals {/j} are given then 
one can find the special contact as A = Xq/Xq{w) where Ao is a contact 
form and w is a contact vector field. Then one can compute the commut- 



ing vector fields v,u\,. . . ,u n -\ from Lemma 3.7 The contact vector field 



defining the dynamical system is a linear combination of these commuting 
fields with constant coefficients. It remains to refer to a theorem of S. Lie 
that if X\ , . . . , X n are commuting and linearly independent vector fields 
in a domain in W 1 then the differential equation x = X\{x) is solvable in 
quadratures, see [3]. □ 
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Remark 3.9 The above theorem is a manifestation of a general phenomenon 
that the existence of an (explicit) Reaction (and even the semi-direct prod- 
uct action) on a manifold implies solvability in quadratures, see O [12] and 
references therein. 

Note that the semi-direct product action, defined via the foliation (J 7 , Q), 
does not depend on a contact form, while the definition of the IR^-action 
requires the knowledge of the special 1-form. 



3.4 Example: integrable geodesic flow on a Riemannian man- 
ifold 

The following is a familiar example from Riemannian geometry. 

Let N be a Riemannian manifold of dimension n, T*N its cotangent 
bundle, H : T*N — ► R the energy function: H(q,p) = \p\ 2 /2, where p is 
the momentum and the norm is the Riemannian one. The geodesic flow on 
T*N is the Hamiltonian vector field of the function H with respect to the 
canonical symplectic structure of the cotangent bundle. 

Note that T*N with the zero section deleted is the symplectization of 
the space M = ST*N of oriented contact elements in N. The homogeneous 
of degree one Hamiltonian \/2H = \p\ defines a contact vector field in the 
contact manifold M; this is the geodesic flow on the space of contact ele- 
ments. The Riemannian metric provides a section of the symplectization 
P = T*N - N -> ST*N = M and hence a contact form on M; namely, M 
is identified with the hypersurface H = 1. The geodesic flow on T*N being 
restricted to M becomes the Reeb vector field of this contact form. 

Assume that the geodesic flow on T*N is completely integrable: there 
exist almost everywhere independent and Poisson commuting homogeneous 
functions /i,...,/ n _i : T*N — > M, invariant under the flow of sgrad H. 
Restricting to the hypersurface M = {H = 1}, one has a co-Legendrian 
foliation J 7 , defined by the integrals and the respective Legendrian foli- 
ation Q, spanned by the Hamiltonian vector fields sgrad /j. The Reeb field 



is tangent to J-, which takes us to the situation of Section 3.2 Thus this 
geodesic flow is a completely integrable continuous time contact dynamical 
system. 

The example of this section can be generalized as follows. Let M be a 
contact manifold, T and Q a co-Legendrian and the Legendrian foliations. 
In the notation of Section |2.2| assume that H : P — > M is a homogeneous 
function of degree one (replacing energy by the norm) , which Poisson com- 
mutes with the functions /». Then the level hypersurface {H = 1} is a 
section of the bundle it : P — > M, and we identify M with this section. 
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Recall that P has the Euler field E, the symplectic structure u and the 
1-form A satisfying relations Q. Let v = sgrad H. 

Lemma 3.10 The vector field v is the Reeb field of the form A. 

Proof. One has: 

i v co = -dH = 0, AO) = u(v, E) = dH{E) = H = 1, 

the first equality due to the fact that if is 1 on the section, and the second 
to the fact that H is homogeneous of degree one. □ 

3.5 Previous work and generalizations 

As we already mentioned, contact complete integrability was studied earlier 
by a number of authors. Here we briefly survey these works. 

P. Liberman [15 a studied Legendrian foliations of contact manifolds en- 
dowed with a contact form A. Such a foliation, Q, is called X- complete if 
the Jacobi bracket of two integrals of Q is again an integral (this does not 
exclude constants). This assumption implies that there exists a flag of foli- 
ations (JF, Q) where T is co-Legendrian and tangent to the Reeb field of the 
form A. It is also proved in [TB] that, in this case, the leaves of T and the 
leaves of Q have affine structures. In our terms, the former is a consequence 



of the local M n -action by contactomorphisms described in Lemma 3.7, and 
the latter is a particular case of Lemma |2.4[ Independently, Pang obtained 
similar results in [19J. 

In terms of the present paper, A. Banyaga and P. Molino [51 El [7] define 
a completely integrable contact manifold as a co-Legendrian foliation whose 
leaves are the orbits of an abelian Lie algebra q of contact vector fields. 
(This point of view was also taken in [IB].) It is proved in [BJ that there 
exists a g-invariant contact form (note that no assumption on compactness 
of the respective group of contactomorphisms is made), and that the Reeb 
field of this contact form belongs to g. Thus one has the situation of Section 



E. Lerman [TH] studied contact toric manifolds, that is, contact mani- 
folds M 2 ™" 1 with an action of a torus T n by contactomorphisms. This is 
analogous to the much better studied theory of symplectic toric manifolds, 
see, e.g., [I]. Let us emphasize that the papers [151 [121 El [BJ El H31 [IB] contain 
many other interesting results; we have mentioned only what is relevant to 
the present work. 
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Remark 3.11 We also note that weakly integrable contact systems dis- 
cussed in Section 2.4 allow various generalizations to manifolds with more 
general non-integrable distributions. Consider a manifold M with a non- 
integrable distribution r, which is not necessarily contact. One way to define 
a dynamical system is to consider on such a manifold the action of a semi- 
direct product group M. 1 ix M fc , which will lead to solvability in quadratures. 
Another way is to consider a foliation trasversal to the distribution r, whose 
intersections with this distribution have a natural M fc -action. This is the 
case, in several examples of non-holonomic mechanics, including the non- 
holonomic oscillator and the Chaplygin skate, which exhibit a weak form of 
integrability, see [5]. 



4 Null lines and the billiard ball map 

4.1 Contact space of oriented light-like lines 

The space of oriented lines M 2n in M n+1 has a canonical symplectic struc- 
ture, which can be defined as follows (see, e.g., |2|). Start with the cotangent 
bundle T*M. n+1 , and consider the unit energy hypersurvace \p\ 2 = 2. The 
restriction of the canonical symplectic structure on T*IR n+1 to this hypersur- 
face has a one-dimensional kernel. The integral curves of this field of kernels 
are called the characteristics. A characteristic consists of unit covectors 
whose foot points belongs to a fixed line and whose kernels are orthogonal 
to this line and agree with its orientation. The space of characteristics is 
again symplectic, and it is identified with the space M of oriented lines in 
R n+1 . 

This construction is called symplectic reduction. Symplectic reduction 
also applies to the space of oriented non-parameterized geodesies of a Rie- 
mannian or Finsler manifold (assuming this space is a smooth manifold, 
which is always the case locally). 

Consider now pseudo-Euclidean space M p,q with p + q = n + 1. There 
are three types of lines: space-like, time-like, and light-like, depending on 
whether the energy \p\ 2 /2 is positive, negative or null. Denote these spaces 
by M? n ,M^" and Mg n_1 , respectively. Symplectic reduction on the energy 
levels ±1 yields symplectic structures on spaces M±, but the symplectic 
reduction on zero energy level yields a space P 2n which is different from 
Mo: the condition \p\ 2 = still allows to multiply p by a non-zero real. 
P is the space of scaled null geodesies which fibers over Mq with fiber M*. 
Thus Mo is a contact manifold whose symplectization is P, the symplectic 
reduction of T*W ,q on zero energy level, see [TO] for details. 
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The space of oriented light-like geodesies was studied about 30 years 
ago by Yu. Manin in his work on application of twistors to the Yang-Mills 
equation. Manin called this space paradise (because it consists of celestial 
spheres, the world lines of photons emanating from point sources in the 
Minkowski space IR 1 ' 3 ). 

4.2 Billiard ball map and accelerating orbits 

The billiard dynamical system in a Riemannian manifold with a smooth 
boundary describes the motion of a free mass-point ("billiard ball"). The 
point moves along a geodesic with constant energy until it hits the boundary 
where the elastic reflection occurs: the normal component of the velocity 
instantaneously changes sign whereas the tangential component remains the 
same. This is the billiard flow, a continuous-time system. The billiard ball 
map T acts on oriented geodesies and takes the incoming trajectory of the 
billiard ball to the outgoing one. T preserves the symplectic structure on 
the space of oriented geodesies. We refer to [201 [21] for information about 
billiards. 

This description applies equally well to billiards in pseudo-Riemannian 
manifolds, in particular, pseudo-Euclidean spaces. A new feature is that 
now the normal vector to the boundary of the billiard table may be tangent 
to the boundary; the billiard reflection is not defined at such points. T 
preserves the type of a billiard trajectory, space-, time-, or light-like. On 
the spaces M±, the billiard ball map is still symplectic, but on the space 
Mo, it is a contact transformation, see |10j . 

In fact, we also have a billiard transformation T : P — > P on the space 
of scaled light-like lines described by the reflection law in the opening para- 
graph of this section. For the projection tt : P — > Mq, one has a commutative 
diagram: tt o T = T o tt. 

Example 4.1 The simplest example is the billiard inside a convex smooth 
closed curve 7 in the Lorentz plane M 1 ' 1 . There are two null directions, 
say, horizontal and vertical, and the billiard system, restricted to the null 
directions, is the following self-map of 7: choose a point x E 7, draw the 
vertical line through x until its second intersections with 7 at point y, draw 
the horizontal line through y until its second intersection with 7 at point 
z, etc., see figure U This map was studied in various contexts, see j9 a for 
references. 

Let us describe an interesting feature of this billiard system, absent in 
the Euclidean case. Suppose we have a closed light-like billiard trajectory. 
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X 



Figure 2: A map of an oval 

Is it possible that, traversing this trajectory, the billiard ball returned to the 
original position with a different velocity vector, say, greater than the origi- 
nal one? Let us call such a conjectural light-like periodic orbit accelerating. 

Assume that the metric is dxdy, so the null directions are vertical and 
horizontal. Let P\, . . . , P2 n £ 7 be the consecutive reflection points of a 
periodic light-like billiard trajectory, and let ti be the slope of the curve 7 at 
point Pi. Consider the billiard ball starting at Pi with, say, unit horizontal 
velocity, (1,0). Then it will return to point P\ with velocity (v,0). 

Lemma 4.2 One has: 

*2*4 • • • hn 
V = . 

tih ■ ■ ■ hn-l 

Further, v = 1 if and only if the periodic light-like trajectory is stable in the 
linear approximation. 

Proof. Consider an instance of reflection, see figure [3] If the slope of 7 at 
the reflection point is t then the tangent vector to 7 is (1, t), and the normal 
vector is (1, — t). Then the reflection is as follows: 

(1,0) = 1(1, t) + \{i-t) ^ \{l,t) - \{i-t) = (0,i). 

Likewise, the vertical-to-horizontal reflection scales the speed down by the 
slope. This implies the first claim of the lemma. 
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Figure 3: Reflection at a point 



For the second claim, consider an infinitesimal horizontal beam reflecting 
in 7, see figure [3] again. Let wi and W2 be the widths of the incoming and the 
outgoing beams. Then, from elementary geometry, W\/w2 = t. It follows 
that the condition for the width of the beam to remain the same after all 
2n reflections is v = 1. But the former is the linear stability condition for 
the periodic orbit. □ 

Since the slopes of 7 at points Pi can be deformed at will (which does not 
affect the reflection, since the orbit is formed by the same null segments), one 
can easily construct a billiard table with an accelerating light-like periodic 
orbit. For such a billiard, there exists no section of the symplectization 
7r : P — > Mo, invariant under the map T : P — » P. 

Remark 4.3 A similar acceleration phenomenon is possible for closed light- 
like geodesies on pseudo-Riemannian manifolds. 



4.3 Billiard inside an ellipsoid 

In this section, we turn to our main example, the billiard system inside an 
ellipsoid in pseudo-Euclidean space. This system was studied in detail in 
|10j : below we summarize the relevant results. 
Consider pseudo-Euclidean space 

yn+l = W ,q w[ih p + g _ n + ^ &nd 

let E : V — > V* be the self-adjoint operator such that the metric is given by 
E(x) ■ x where dot denotes the pairing between vectors and covectors. Let 
A : V — > V* be a positive-definite self-adjoint operator defining an ellipsoid 
A{x) ■ x = 1. Since A is positive-definite, both forms can be simultane- 
ously reduced to principle axes, and we assume that A = diag(a^ 2 , . . . , a~ 2 ) 
and E = diag(l, . . . , 1, —1, . . . , —1). Consider the pseudo-confocal family of 



20 



quadrics 



^ + ^ + --- + ^ + ^ ±L T + --- + ^ X T = 1 (3) 
af + A a\ + A a£ + A a^ +1 - A a^ +g - A 

where A is a real parameter (see figure [4] for a two-dimensional example) . 
Let Mq™" 1 be the contact space of oriented null lines in V n+1 , and let P 2n 
be its symplectization, the space of scaled null lines. 



Figure 4: A family of pseudo-confocal conies; null directions have slopes ±1 

The following theorem is proved in |10j . 

Theorem 4.4 1 ) The tangent lines to a fixed light-like geodesic on an ellip- 
soid in pseudo-Euclidean space V n+1 are tangent ton— 2 other fixed quadrics 
from the pseudo-confocal family |3p. 

2) A light-like billiard trajectory inside an ellipsoid in pseudo- Euclidean 
space V n+1 remains tangent to n — 1 fixed pseudo-confocal quadrics. 

3) The set N of oriented light-like lines, tangent to fixed n—\ pseudo-confocal 
quadrics, is a codimension n — 1 submanifold in Mq, foliated by Legendrian 
in Mq submanifolds, which are of codimension one in N. 

(For space- and time-like lines, the number of pseudo-confocal quadrics in 
statements 1 and 2 is one greater.) 



21 



In terms of the present paper, the null lines, tangent to n — 1 fixed 
pseudo-confocal quadrics, constitute the leaves of a co-Legendrian foliation 
T in Mq. Let /i,...,/ n _i : Mq — > R be smooth functions defining the 
foliation .F (these functions index the pseudo-confocal quadrics tangent to a 
given line), and let fi be their lifts to P, the space of scaled light-like lines. 
Then the functions fi Poisson commute. 

Let us also describe the leaves of the Legendrian foliation Q. A leaf of T 
consists of null lines I tangent to fixed n — 1 pseudo-confocal quadrics, say, 
Qi, ■ ■ ■ , Qn-i- Let Vi be the geodesic vector field on TQi. Considering the 
oriented tangent lines to a geodesic curve, we view as a vector field on the 
space of lines tangent to Qi. Then these vector fields commute, and the leaf 
of the foliation Q through point I is generated by the fields v±, . . . , v n -\. 

Explicit formulas for integrals are as follows, cf. |17|[T%] in the Euclidean 
case. Identify the tangent TV and cotangent T*V spaces via the pseudo- 
Euclidean metric. Then one has the following integrals of the billiard flow 



where the coordinates of the position and Vi of the velocity vectors, 

and where e\ = ■ ■ ■ = e p = 1, e p+ \ = • • • = e p+q = —1. These integrals 
Poisson commute and satisfy the relation ^ = (v, v). The same formulas 
give integrals of the geodesic flow on a quadric in pseudo-Euclidean space 
V. Note that the integrals Q are quadratic in velocities. 

In the Euclidean case, when all e< = 1, the functions Fk/(v,v) descend 
to the space of oriented lines and are integrals of the billiard ball map. In 
the pseudo-Euclidean case, (v, v) = for the null directions, and one cannot 
divide by (v, v). 

Following |21j . let us describe another integral of the billiard ball map, 
homogeneous of degree one in the velocity. Let x be a point of the ellipsoid 
and v an inward vector with foot point x. As before, one has the billiard 
ball transformation T on such tangent vectors. If v is null then the set of 
the inward tangent vectors with foot point on the ellipsoid is identified with 
the space of scaled oriented lines P. 

Proposition 4.5 1) The function H(x, v) := Ax ■ v is negative. 

2) H(x,v) is invariant under the billiard ball transformation T. 

3) H(x, v) Poisson commutes with the functions fi, i = 1, . . . ,n — 1. 

Proof. For the first claim, note that Ax is the outward normal covector 
and v has the inward direction, hence H (x, v) < 0. 



on TV: 
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For the second claim, the billiard ball map is the composition of two 
maps: (x,v) \— * (y,v) i— > (y,u), where the second is the billiard reflection, 
see figure [5] We claim that Ax ■ v = —Ay ■ v = Ay ■ u. 





Figure 5: Billiard ball map as the composition of two involutions 

To prove the first equality, note that {Ax + Ay) -{y — x) = since A* = A 
and Ax ■ x = Ay ■ y = 1. On the other hand, v is collinear with y — x, hence 
Ax ■ v = —Ay ■ v. 

To prove the second equality, note that, due to the reflection law, v + u 
is a tangent vector to the ellipsoid at point y. On the other hand, Ay is 
the normal covector to the ellipsoid Ay ■ y = 1. Hence Ay ■ v = —Ay ■ u, as 
claimed. 

Now we prove the third claim. Extend the function H to the tangent 
bundle TV in such a way that it is invariant along straight lines: H(x + 
tv, v) = H(x, v) for all t £ M. Then this extended function is an integral 
of the billiard flow inside the ellipsoid since it is also invariant under the 
reflection, see above. Since is a complete system of first integrals, H 
is functionally dependent on integrals in Q. Note that each is also 
invariant along straight lines: Fj-(x + tv, v) = Fj-(x,v). Hence the functional 
relation descends, in particular, to the space of scaled null lines P. 

Therefore, it suffices to show that the functions ft and F^ Poisson com- 
mute in the space P of scaled light-like lines. Indeed, as we mentioned 
above, the Hamiltonian vector field Vi = sgrad/j defines the geodesic flow 
on the pseudo-confocal quadric Qi, and the functions F^ are integrals of the 
geodesic flow on these quadrics, hence {/j,-Pfc} = 0. □ 



Proposition 4.5 places us in the situation of Section 3.4 and therefore, of 



Section 3.3: one has an invariant contact form on the space of null geodesies 



whose Reeb field is tangent to the co-Legendrian foliation. More specifically, 
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the integral H (x, v) provides a section of the symplectization bundle P — > 
Mo, and hence a special contact form on Mq. It follows that the billiard ball 
map on the light-like oriented lines inside an ellipsoid in pseudo-Euclidean 
space is a completely integrable contact transformation. 

Remark 4.6 One can view the contact integrable system on null geodesies 
on an ellipsoid as a limit of the Hamiltonian completely integrable system of 
space-like geodesies. While only n — 1 independent integrals of the Hamilto- 
nian system (out of the n) survive in the limit, when passing to the contact 
manifold, the corresponding R"-action on the symplectic manifold of space- 
like geodesies does extend to the space of null geodesies on the ellipsoid. 
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